Strength functions, entropies and duality in weakly to strongly interacting fermionic 

systems 



I. INTRODUCTION 



D. Angom, S. Ghosh, and V.K.B. Kota * 
Physical Research Laboratory, Ahmedabad 380 009, India 

' We revisit statistical wavefunction properties of finite systems of interacting fermions in the 

' light of strength functions and their participation ratio and information entropy. For weakly 

, interacting fermions in a mean-field with random two-body interactions of increasing strength A, 

(N , the strength functions Fk{E) are well known to change, in the regime where level fluctuations 

(-H ' follow Wigner's surmise, from Breit-Wigner to Gaussian form. We propose an ansatz for the 

, function describing this transition which we use to investigate the participation ratio ^2 and the 

information entropy 5""'° during this crossover, thereby extending the known behavior valid in 
' the Gaussian domain into much of the Breit-Wigner domain. Our method also allows us to derive 

the scaling law for the duality point X — Xd, where Fk{E), £^2 and 5""^° in both the weak (A = 0) 
and strong mixing (A = 00) basis coincide as Ad ~ l/v^m, where m is the number of fermions. As 
an application, the ansatz function for strength functions is used in describing the Breit-Wigner 
^ ' to Gaussian transition seen in neutral atoms Cel to SmI with valence electrons changing from 4 to 8. 
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There are many physical systems which are statistically well described by so-called embedded ensembles of fermions, 
C I representing particles subjected to a one-body mean-field potential (defining a set of single-particle levels), and 

interacting with a random, two-body potential. Examples include heavy nuclei Q, |^ H, Q , natural 0, , or artificial 
atoms (quantum dots) 0, I3j ^ind nanometer-scale metallic grains Similar situations of randomly interacting 
spin systems occur in the study of spin-glass systems and in the context of quantum information and quantum 
computation These embedded ensemble are defined as ensembles of Hamiltonians {H} = h{l) + X{V{2)}, 

where {. . .} denotes an ensemble distribution, /i(l) = eiUi is a fixed one-body operator (one can also consider an 
;J| ' ensemble {/i(l)} defined by a probability distribution P{ei)) defined by the single-particle energies with average 

5^ ! spacing A which sets the energy scale (one can thus set A = 1 without loss of generality), and n.i is the number 

operator for the single-particle state \i). Similarly V{2) is the random two-body interaction with its two-particle 
matrix elements chosen as independent Gaussian variables with zero center and unit variance. Thus, for m fermions 
in TV single-particle states, {H} is a oneplus two-body random matrix ensemble (called embedded ensemble of 
(l+2)-body interactions[EGOE(l-|-2)]) PiQl defined by the parameters (m, N, A), where A is the interaction strength. 
For convenience, we only consider here EG0E(l-|-2) f or sp inless fermions, but extensions to particles with intrinsic 
angular momentum have also been considered 0, M In such a case, the size of Hilbert space is given by 

^ im) another important parameter is the connectivity K giving the number of directly coupled states as 
A' = 1 -I- m{N -m) + m{m - 1){N - m){N - m - l)/4. 

Because of its broad relevance to many, a priori different, finite quantu m sy stems, EG0E(l+2Vs have been 
investigated in detail by many research groups in the recent past [1 11 |lB II II 111 111 El 111 E El 111 111 111 113 ■ 
Most investigations used analytical methods extrapolating from the weak and the strong coupling limit, and relied on 
numerical calculations in the regime of intermediate values of A. Focusing on the statistical spectral and wavefunction 
properties, the dominant features of EG0E(l+2) that emerged from those investigations can be summarized as: 

1. There is a marker Ac, such that for A > Ac the many-body level spacing distribution becomes close to that of the 
Gaussian Orthogonal Ensemble (GOE) of random matric es l2lll. while for A < Ac the level fluctuations are close 
to Poisson. It was further established that Ac oc l/m?N [lajlla; specifically for m = 6 and N — 12, Ac — 0.06 

Q. 

2. As A increases from A = 0, the strength function (to be defined below in Section 2) undergoes a crossover 
from a delta-peak, first to a Breit-Wigner (BW), then to a Gaussian form. Related to that crossover, there are 



2 



two markers A^^'^-* such that, for A^^"* < A < Xp\ the strength function is well approximated by a BW form 

P, U IT^ IT^ ITtI [T9I I ; the BW form emerges above Ay^'' , which is exponentially smaller in m and N than 

Ac- In particular, the BW form occurs while the spectral fluctuations are still Poissonian. The X > Xp region 
(with full GOE spectral and wavefunction properties) is called the Gaussian domain Q . From now on we put 

A_F = A^^; note that Xp » Ac. Arguments based on BW spreading widths give Xf oc [mil; for the 

specific case m — 6 and = 12, Af — 0.2 Q. 

3. In the Gaussian domain, the participation ratio (PR) ^2 and the exponential of the information entropy 
(exp[S''"^°(i?)]) (both quantities will be defined in Section 3 below) take Gaussian forms when plotted as a 
function of the energy Q. The variances of these Gaussian are (1 + C^)/(2C^) and 1/C^ respectively, where 

= f^(w)/[tT^(TO) + A^cTy(m)]; (Jh{m) is the spectrum width produced by h{l) in the total m-particle space 
and similarly cry (to) is the width produced by 1^(2). Also, in the BW region, the PR is given by the ratio of the 
spreading width and the spacing between directly (by V{2)) connected states and S™^° ~ In(PR) (see also 

4. There is a fourth marker Ad such that at A = Ad the strength functions, PR and 5*'"^° expressed in either the 
h{l) (i.e. A = 0) and V{2) (i.e. A = 00) basis will coincide. This is accompanied by a duality transformation 
relating the values of those quantities in the h{l) basis to those in the V{2) basis by A ^ X'^/X 0. In Section 
5 the (to, A') dep endence of Ad will be shown to be Ad ~ iZ-^ro (correcting the previously postulated result 
Ad ~ 0); for m = 6 and A^ = 12, as we see ahead. Ad — 0.3. 

In addition, another important result for EG0E(l+2) is that the smoothed (ensemble averaged) density of states takes 
on a Gaussian form independently of A [lll^l23l|. Our purpose in this paper is to bring completion to the investigations 
related to the points (l)-(4) above. In particular, we will introduce an interpolating function for strength functions 
for the BW to Gaussian transition and apply it to neutral Cel to SmI atoms. Now we will give a preview. 

In Section 2 we discuss a variant of the well known Student's t-distribution (hereafter called Fk:BW-g{E)) and 
show that it is well suited for describing the BW to Gaussian transition, as a function of one parameter a. Numerical 
calculations allow to establish a one-to-one correspondence between a and the interaction strength A. In Section 3, 
the resulting Fi^.Bw-g{E) is used to calculate both PR and 5'"*°, and comparison is made with direct numerical 
calculations of these quantities as a function of A, over the full range of variation of A, thereby extending previous 
similar investigations which were restricted to either the Gaussian ^ or BW 0,^3 domains. Additional structures 
in the wave functions can be captured by the structural entropy Sst-c = S™^° — In ^2, which measures the amount of 
information contained in the tails of the strength functions. Results of an analysis of Sst-c is also given in Section 3. 
In Section 4, the Fk:BW-g{E) is applied in the analysis of the BW to Gaussian transition one observes as we go from 
neutral Gel atom to SmI atom. In Section 5, the existence of a duality transformation in EG0E(l-|-2) (which was 
established in Ref. ^1) is discussed and it is shown that, using the results for PR and 5'"^° in the Gaussian domain, 
the duality point Ad ~ \/\/rn. Conclusions and final comments are given in Section 6. 



II. INTERPOLATING FUNCTION FOR BW TO GAUSSIAN TRANSITION IN STRENGTH 

FUNCTIONS 

Given the mean-field h{l) basis states |fc) and the expansion of the eigenstates \E) as \E) = \k), the 

strength function Fk (E) is defined by 

F,(i?) = 5^ |cf |'<5(i?-i?')-(|Cfr> {dp"{E)) (1) 

E' 

In Eq. (1), (...) indicates an ensemble average, d = (^) is the TO-particle space dimension and [E) is the 
normalized (and ensemble averaged) density of states. As mentioned in Section 1, p^{E) is in general a Gaussian 
(often the superscript H is dropped), 

p"iE)= ^\ ^ exp-:^; E = (E - eHim))/aHim) (2) 
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where euirn) — (H)"^ is the spectrum centroid and similarly (Ti/(m) is the spectral width. The BW and Gaussian 
(denoted by Q) forms of F^iE) are, 

where Ek — {k\H\k). With p = J^^^ Fk{E)dE, the spreading width — £3/4 — £1/4. Similarly the variance of 

Fk is a'j, — (k\ \ fc) — ((fc | H \ k))^. Expressions for Tk and ak in terms of (m, iV, A) are given in Sections 5 
and 3 respectively. Both spreading width of the BW and ak of the Gaussian strength functions are essent ially 
independent of k (provided one considers energies not too far away from the center of the density of states) 0, 0, H^ . 
Simultaneously, the energies E^s (of H) and EkS will have the same centroid. Moreover, just as the state density 
p^{E), the EkS density (denoted by p^^[Ek)) is also a Gaussian. These results are used throughout this paper and 
without loss of generality the centroids of E's and EkS are set equal to zero. As it is discussed in detail in Q, the 
width of p^^{Ek) is essentially due to h{l) (with a small correction from V{2) as explained in the Appendix in 
and (Jk is generated by V{2). Before proceeding further it should be mentioned that the strength functions are basis 
dependent and one can define strength functions in the V{2) basis also. We will turn to this question when discussing 
duality transformation in Section 5. 

For the BW to Gaussian transition we make the following ansatz for Fk{E) 

Fk..BW-Q{E:a,P)dE=^ Y , a>l 4 

The Fk:BW-g in (4) gives BW for a = 1 and Gaussian for a 00 (this can easily be checked using Stirling's 
approximation). As required, it is normalized to unity for any positive value of the continuous parameter a. For 
2a — 1 a integer, Fk-BW-g gives the so called Student's ^-distribution |24| . which is well known in statistics. In 
particular, the Student's distribution f{x) with a parameter v given in Table 5.7 of reduces to (4) with the change 

a — {v + l)/2, V a positive integer, and x ^ {E — Ek)/y/p. Note that the construction of Fk-.BW-g in Eq. (4) 
is similar in spirit to the Brody distribution, interpolating between the Poisson and Wigner-Dyson distributions for 
nearest neighbor spacing distribution (NNSD) 2]. Also, just as some groups use for the NNSD a linear combination of 
Poisson and Wigner forms multiplied by x and (l — x) respectively with x being the mixing parameter, it is possible to 
use pFk:Bw{E) + (1 — (x)Fk:g{E) for the BW to Gaussian transition with /i (0 < /x < 1) being the mixing parameter. 
This simple form is not explored in this paper as it is unlikely that a theory for strength functions for EG0E(l+2) 
will give this form. 

In Fk;BW-g{.E : a, (3), the parameter a is sensitive to shape changes, while the parameter /? supplies the energy scale 
over which Fk:BW-g{E : a,/?) extends. Since we focus on the shape transformation, a is the significant parameter. 
First, it is easy to see that Fk-BW-g{E : a, (3) is an even function of E — Ek, so that all of its finite odd cumulants 
vanish (strictly speaking, the centroid is Ek only for a > 1; see 24] ). The variance of Fk.BW-g, defined only for 
a > 3/2, is 



cT'{Fk..Bw-g) = [2^^) " > 3/2. (5) 

and it is useful to recall that a'^{Fk) — X^<Jy For a > 3/2 one can use (5) to fix /3 while for a < 3/2, it is 
the spreading width F (this is well defined for all a values) that is useful for fixing the f3 value. There is no simple 
expression for F as a function of a and f3 but using (4) this can be calculated numerically. Just as Eq. (5) , the excess 
parameter (also known as Kurtosis) of Fk.BW-g is 72 — 6/ {2a ~ 5) for a > 5/2. However this expression is not useful 
in practice as the spectrum is always of finite range and this causes large deviations for a « 2 — 8. Therefore it is 
more useful to consider 72 with the spectrum ranging say from —a to +a. Then Eq. (4), with proper normalization, 
gives, 

, 9 2fi(|,a;i;-77^)2fi(ia;|;-772) 

72(a :«,/?) = - \ ; ^-72^ 3 (6) 

5 2i^i(f,a;f;-?72) 



4 



where rf — {a'^ /a/3), and 2F1 a hypergeometric function. 

In Fig. 1 the results of EG0E(l+2) for Fk{E) are compared, for the m = 6, iV = 12 system at — (i.e. in the 
middle of the band), with the best fit Fk-.BW-g for various values of A. In the fits, the f3 values are fixed using Eq.(5) 
for A > 0.1 (for these, a > 1.6) and the spreading width F is used {(3 ^ F'^/4) for A — 0.06 (here a = 1.2). In the fits, 
also imposed is the condition that the value of 72 calculated from Eq. (6) over the spectrum range should be close to 
the numerical EG0E(l+2) values. As it is seen from Fig. 1, the fits are excellent over a wide range of A values; in the 
calculations only A > 0.06 are considered (for the system considered in Fig. 1, Ac ~ 0.06). The variation of a with 
A is shown in Fig. 2. The parameter a raises slowly upto A^ (note that Xp ^ 0.2 for the EG0E(l+2) system used 
in Fig. 1 01) and then it starts rising sharply with A. Finally the a values start saturating after A > Aq = 0.3 (the 
saturation is artificial as determination of a for A > > Aq is difficult as here Fk {E) will be very close to Gaussian) . The 
criteria a ~ 4 and 72 '--^ 1 appear to define Xp- Fig. 2 shows that the BW to Gaussian transition is a sharp transition 
and therefore studies in BW and Gaussian regimes can be carried out independently, to a good approximation, as it 
is done in many papers before. Now we will apply Fk-Bw-g to study PR and S*™-^" in the region intermediate to BW 
and Gaussian forms. 



III. PARTICIPATION RATIO AND INFORMATION ENTROPY IN THE BW TO GAUSSIAN 

TRANSITION REGION 

Two important measures of the complexity of eigenstates of interacting systems are the participation ratio and 
information entropy. As in the previous section, we expand the Hamiltonian eigenstates in the noninteracting mean- 
field basis as \E) — \k). Then PR and information entropy are defined as 

-1 

(7) 



The subscript '2' denotes that ^2 is the second Renyi entropy [23 • Qualitatively, ^2 counts the number of {|fc)}-basis 
components necessary to construct one typical |i?)-state, and is thus often referred to as the Number of Principal 
Components (NPC). Obviously, both ^2 and S™^° are basis dependent, and could as well be defined starting from 
another expansion. Eq. (7) gives their expression with respect to the h{l) basis (i.e. the |fc)'s in Eq. (7), as in 
Section 2, are h{l) basis states), consequently, ^2 and give measures of the spreading of eigenstates over the 

noninteracting basis as the many-body interaction is made stronger and stronger. Below, in Section 5, we will deal 
with these measures defined with respect to the V{2) basis. Following Ref. ,J| one can write ^2 and S™^° in terms of 
the strength functions. 




{UE)/i§^''y' = ——2 / dEkp"^{Ek)[Fk{E)f 



1 f°° 

S'-'°{E)^S'SoE = -^TTin / dE^ p"HEk)FkiE) In 



Fk{E) 
p"{E) 



(8) 



An obvious notation refers to the GOE values. 



(9) 



In the Gaussian domain, with Fk{E) being a Gaussian with centroid at E^ and width XaviiTi), the integrals in Eq. 
(8) are easy to evaluate (note that p^{E) and p^^{Ek) are zero centered Gaussians with width a Him) and ah{m) 
respectively) and they give Q, 



Ce^ 



6(i?)/e2°°^ = v/W^cxp- 



exp(^-f°(£;) - ^j^g'fi) = VW^exp^ exp 



(10) 

C Ce^ 



2^2 
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The important parameter in Eq. (10) is the correlation coefScient (. For ^2 and 5""^° in h{l) (A = 0) basis C = Co 
where 



(Am)\ ' ^ ^'M nil 



Strictly speaking ct^ (m) in the numerator should be replaced by the variance of energies and similarly the denom- 
inator by ajj{m). Eq.(ll) is obtained by recognizing that the former is very close to af-^ and the later is essentially 
cr^ + X^cTy. In fact these results are valid in the dilute limit (m ^ 00, iV ^ 00, m/N 0) and here h and V are 
orthogonal. Even away from the dilute limit they remain to be good approximations (see Fig. 3 ahead for a test). 
Propagation formulas 0] for cr^ (to) and ay (to) are 

(12) 

_ m{m - l)iN - m){N - to - l)jV(iV - 1) _ 
^ (iV-2)(iV-3) ^ 

It is always possible to write c^(l) in terms of and for example for a uniform single particle spectrum, 

alil) = iN+l){N-l)^ (13) 

As shown in Fig. 3 (for the m — 6, N — 12 system), results of the formulas (11,12) agree very well with numerical 
EG0E(l+2) values for Co- 
Substituting the interpolating Fk:BW-g{E) for Fk{E) in Eq. (8), one can study ^2 and 5'"^° as a function of A. 
The integral in Eq. (8) for 1^2 can be simplified for £ = and this gives (for a > 3/2), 



THa) 1 ^A3_ (2«-3)(l 




(14) 



(2a -3) r2(a-i) VC^(l-C^) V2'2 ' 2^ 

where U{ ) is hypergeometric-U function 26]. The corresponding result in the Gaussian domain (from Eq. (10)) 

IS yr^. In der iving Eq. (14), we used Eqs. (5) and (11) for eliminating /3 and simplifying all the variances 
into C^. It is important to note that Eq. (14) is valid only for a > 3/2 as (5) is used in deriving this formula. For 
a < 3/2 a compact formula could not be derived but one can use (8) for numerical evaluations. Similarly, in the case 
of S"^^°{E = 0) a simple formula like Eq. (14) could not be obtained for any a and once again here one can use (8) for 
numerical evaluations. For £,2{E — 0) results from (14) for A > 0.08 and the result from (8) for A = 0.06 are compared 
with numerical EG0E(l-|-2) calculations for the m = 6, 12 system in Fig. 4. In these calculations the a-values 
are read off from Fig. 2 and C^ from Fig. 3. Comparing with the Gaussian domain results given by (10), it is seen that 
they are good for A > Af as expected; these results again confirm that Xp ^ 0.2 for the m = 6, A*" = 12 system. The 
agreement between (14,8) and the numerical calculations continue upto A ^ 0.06. For A < Ai? as the BW structure is 
more dominant, there will be more localization and hence ^2 decreases fast as A is decreasing as seen in Fig. 4; finally 
they will approach zero for A ^ 0. The results based on (8) will not extend to the region A ^ Ac as here the GOE 
assumptions used in deriving these equations (see Q) will fail. Finally, for S™^°{E — 0) the results are similar to those 
shown in Fig. 4. This is not surprising as in many numerical calculations (including the present calculations) it is seen 
that S"^^°{E) ~ ln(^2(£')) and therefore only their difference can capture the information not contained in the bulk 
of S""f° or PR. With this clue, recently it is argued [H that the structural entropy Sstr{E) = S'"^^°{E) - ln{£,2{Ej) 
is an important measure of complexity (in addition to S'™^°(i?) or £,2{E)) in eigenfunctions. More importantly S'str is 
free of divergences associated with S*'"*^" and PR. Note that exp(S'QQ^) and £,2^^ are 0.48c? and d/3 respectively and 
therefore they diverge as the matrix dimension d —^ 00. For interacting particle systems it is observed that S'str (-E = 0) 
vs A (or the disorder in the Anderson model [25l|') exhibits a peak. It is then of interest to examine S'str in terms of the 
results given Section 2. First we consider Sstr in the BW and Gaussian domains. 
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For small A, one can estimate Sstr using the the BW approach, i.e. taking Fk-^w from Eq.(3), and using (|C|'p) = 
Fk:Bw{E^)^rrn where gives the many-body level spacing. Inserting this into Eq.(7) and replacing the sums by 
integrals over one then gets, by restricting to £' = 0, 

- 1^ , a = arctan(2B/r) (15) 

where B is the many-body bandwidth and F (x is the BW width. Fhe upper limit 5str(0) ^ ln2 0.7 follows by 
letting B/r oo (then a 7r/2) in (15). With F increasing with increasing A, the 5811(0) starts decreasing from the 
maximum value. Similarly, in the Gaussian domain, using (10), one has 

S^tM = ln(1.44) + ln(l + C')) (16) 

It should be noted that S^^^ ~ ln(1.44) independent of E. An interesting observation (though its significance is 
not clear) is that for A = (then C = 1) the Ssti is sum S'str for GOE and a Gaussian; note that for a Gaussian 
distribution, as shown in (2^, 5str = ^(1 — ln2). As A ^ cx) gives C = 0, S'str starts from GOE value (0.3689) for very 
large A and then starts increasing with decreasing A. 

For the intermediate regime and for A < Ac no analytical results could be derived yet but numerical calculations 
give some insight. Fig. 5 shows the EG0E(l-|-2) results for SstriE = 0) vs A in the /i(l) basis and their comparison 
with the results from (8,14) where Fk-.BW-G given by (4) is used. The S'str is well described for A >> Ac and for 
the m = 6, = 12 system considered in Fig. 5 the theory given by (8,14) is good upto A — 0.1. Comparing with 
Eqs. (15,16), it is seen that the Gaussian domain result (16) describes the results for A > Af while the BW result 
(15) describes only the trends for A between 0.1 and Xp (with the maximum possible value for S'str being 0.7). More 
importantly, as seen from Fig. 5 and also from Fig. 4 of |l9j| . Sstr exhibits a peak around a A value not far from 
Ac marker and here the level fluctuations will have Poisson component. Thus it is plausible that the peak arises due 
to large spectral(and strength) fluctuations. A good theory for S'str generating the observed peak is at present not 
available. 



Sstr(O) 



tan a j^.^ 

de- 

' tan a 



1 



In 



2a -t- sin 2a 



In 
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IV. BW TO GAUSSIAN TRANSITION IN NEUTRAL ATOMS CEI TO SMI 



Realistic examples for the BW to Gaussian transition in strength functions are expected to come from neutral 
lanthanide atoms. It is known from the analysis of Cel by Flambaum et al and PrI by Cummings et al that 
the Fk{E) for these atoms with 4 and 5 valance electrons respectively are close to BW while those of SmI with 8 
valance electrons, as shown in are close to Gaussian. This series of atoms is completed by Ndl and Pml with 6 
and 7 valance electrons. We made calculations for not only Ndl and Pml but also for the other three by using the 
same method and this is briefly discussed below before giving the results. Also let us add that as Af cx: l/^/m, it is 
to be expected that, with m changing from 4-8, there should be BW to Gaussian transition with H fixed. This is 
indeed confirmed by the results discussed ahead. 

The ground state configurations of Ce, Pr, Nd, Pm and Sm are Af5d6s'^{^Gi), 4:p6s^{'^lQ/2), if^es^i^h), 
4/^6s^(^i/5/2) and 4/^6s'^(''Fo) respectively. Couphng of the 5c? and 4/ valence electrons produce several config- 
urations and indicate strong configuration mixing. Previous work on Sm I and lanthanide series [29j | have 
established that an appropriate method of calculation is the multiconfiguration Dirac-Fock (MCDF) [s^l, where an 
atomic state function (ASF) \rPJM) is approximated as a linear combination of configuration state functions ( CSFs) 
I7PJM). That is \TPJM) = J^k'^kllrPJM), where c/jS are the mixing coefficients, P, J and M are parity, total 
angular momentum and magnetic quantum numbers respectively, and F and jk are additional quantum numbers to 
define each of the ASFs and CSFs uniquely. The CSFs are linear combinations of Slater determinants and ASFs 
are eigenfunctions of the Dirac-Coulomb Hamiltonian H^'^ [s^l • In the present study, a series of extended optimized 
level(EOL)-MCDF calculations are carried out using GRASP92 to generate a single electron basis set consisting 
of (1 - 6)si/2, (2 - 6)pi/2, (2 - 6)p3/2, (3 - 5)^3/2, 4/5/2 and 4/7/2 orbitals. The EOL-MCDF orbitals are less state 
specific compared to ground state extrcmization and suitable for studying spectral properties and structure of ASFs. 
Then a CSF space of a specific J having single and double excitations from a reference configuration 4/'5d™6s^ to 
5d, 6p and 4/ shells, where I and m are the occupancy of the shells is generated. This sequence of calculations is 
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TABLE L Details of the angular momentum, the reference configuration considered and the number of CSFs generated for 
each of the atoms. The numbers within parenthesis are the number of CSFs chosen for the final calculations. 



Element Ref Config Number of CSFs 

Ce 4- 4f5d6s^ 373 (308) 

Pr 9/2" 4f6s^ 1378 (1278) 

Nd 4+ 4/'*6s^ 2200 (2000) 

Pm 9/2" 4f6s'^ 4378 (4178) 

Sm 4+ 4f6s^ 7325 (6500) 



repeated for each of the atoms. Details of the and the reference configuration considered and the number of CSFs 
generated for each of them are given in Table HI Note that the parity (P) is chosen to be same as that of the the 
ground state and J is 4 for even and 9/2 for odd cases. ASFs and corresponding eigenvalues within the CSF space 
considered are obtained by a configuration interaction calculation. For further analysis related to the structure of the 
ASFs, we choose CSFs which have close to uniform separation of H^jp = {'~fkPJM\H^'~^\jkPJM). For example, in 
Sm only 6500 of the 7325 CSFs generated are considered, the first 200 and last 625 CSFs are excluded. Then, the 
strength function Fk{E) = J2e' I'-fe \'^S{E — E') of the selected CSFs are calculated, where is the mixing coefficient 
of \^kPJM) for the ASF having eigenvalue E. To get a representative Fk{E) we calculate the average of Fk{E)s 
around the centroid and over 3% of the range of Hf^, i.e. average strength functions with Ek — 0. 

Calculated average strength functions (with Ek = 0) are compared with the best fit Fk-BW-g{E) in Fig. 6. In 
the fits, Eq. (5) is used to eliminate /3. The a values for each atom are given in the figure. Firstly it is seen that 
Fk:BW~g{E) gives excellent description of the calculated strength functions and the BW to Gaussian transition is 
clearly seen in Fig. 6 with a changing from 1.85 to 14 as we go from Cel to SmI. The calculated 72 values are also 
consistent with this transition as they change from 6.44 to 0.46. Comparing with Fig. 1, Cel and PrI atoms are 
close to A ~ 0.1 — 0.15 cases {Fk is close to BW), Ndl and Pml are close to A ~ 0.2 — 0.25 cases {Fk is intermediate 
to BW and Gaussian) and SmI is close to A ~ 0.3 case {Fk is close to Gaussian) of the m = Q,N — 12 EGOE 
example. For further confirming the BW to Gaussian transition the ^2(^' = 0)/^2''^^ values are calculated using Eq. 
(14) and the deduced a values (used also are the calculated C values). They change from 0.21 to 0.6 for Cel to SmI. 
These EG0E(H-2) numbers are close to ^2{E — 0)/^^*^^'s generated by the calculated atomic eigenstates. However 
there are large fluctuations in ^2(^)1 as the atomic calculations produce in general, for many states, more localization 
than expected from EG0E(l+2); if we average ^2{E) in the neighborhood of £^ = 0, then the calculated values are 
~ 20 — 30% smaller than the values given by Eq. (14). This is already seen in Cel and PrI in [23| and SmI in |^. 
The present calculations confirm this to be a generic behavior. The source of this localization and modifications of 
EG0E(l+2) for incorporating this property will be discussed elsewhere. Here it suffices to conclude that, from the 
results in Fig. 6, Cel to SmI unmistakingly exhibit BW to Gaussian transition for Ek = 0. 

V. DUALITY BETWEEN WEAK AND STRONG MIXING LIMITS 

For each realization of the Hamiltonian H — h(l) + XV(2), two asymptotic natural basis emerge: the noninteracting 
basis defined by h, and the A = c» basis defined by V. The previous discussions in Sections 2 and 3 were concerned 
with strength functions, PR and 5™^° in the noninteracting basis only. Here we extend this discussion to the A = 00 
basis, and will focus on the existence of a duality transformation between the two basis, following the recent work of 
Jacquod and Varga ( JV) . JV found that a duality point A^ exists where all the statistical wavefunction properties 
in these two basis coincide, and that the wavefunction properties in the noninteracting (A — 0) basis are related to 
those in the X = 00 (fully interacting) basis by a duality transformation A <-> A^/A. An ambiguity in JV results lied 
with the fact that the existence and scaling of the duality point A^ were derived within the BW approximation, while 
Xd explicitly lies outside the BW regime. We therefore extend those theoretical arguments by similar ones in the 
Gaussian approximation, but first recall JV results. 

In the noninteracting basis, and in the BW regime, the strength functions will have Lorentzian shape. To estimate 
its width via the golden rule, one first has to realize that there are, beside the one-body spacing A, two important 
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energy scales ■ the mean spacing between states directly coupled by the two-body interaction Ac"' — 8^2''' / K w 

4A /Nm? and the m-body spacing Am' — Bjn /d, where i?m' ~ ^/mNA is the m-body band (note that Bm^ ~ ah {m) 
with ah{m) given by (12,13)). Note that this estimate slightly differs from that of JV where the m-body bandwidth 
Sm' was approximated by the extremal possible energy values, instead of the rms of the density of states. Then, the 
width of the BW strength function is approximated via the golden rule as r*^"^ cx A^/Ai^' ^ X'^Nrn? / A. Finally, in 
the dilute hmit m < N, the PR is obtained as ^^°' = r("V^™^ « A^m^/^d/A^. Note that, this resuh differs from 
the JV estimate given in (T^ by the factor mP^^ instead of m. As is the case for r(0), the Golden rule gives a good 
estimate of the width r^°°> of the strength functions expressed in the A = 00 basis, in the BW regime. Following JV, 
it is seen that r(°°) - m{N - rn)A'^/X, and one gets the PR in the X = 00 basis as f^°°' = r(°°)/A^' cx {A/Xfd. 
The duality point is then defined by 'C2*''(-Arf) = ^2^^ {^d) which gives the parametric dependence 

Ad cx A/m^/s. (17) 

This result |32j | is in better agreement with the numerical data presented by JV, where it was found that A^ ^ 1/ rri'^ 
with v S [0.3,0.5] (one has to keep in mind however, that most data were not in the dilute limit and that v was 
extracted from a restricted range of variation of m); compared to the previous JV estimate A^ cx A/m}/^, Ea. (|17() is 
thus in better agreement with numerical data. 

A very close estimate for A^ can be derived from PR and S™^° in the Gaussian domain. Just as C = Co in Eqs. 
(10) describes S™^° and PR in the Gaussian domain, 

Co(A) - cjj^al + X^al - ^{PA^)I{PA^ + g^X^) , (18) 

it is to be expected (by extending in a straight forward manner the arguments in j3| where h{l) basis is considered) 
that in the V{2) basis also the PR and 5'"*° will be given by (10) but with C = Coo where. 



C^(A) = Xavl^al + X^al - ^ {g^ X^) / [p A^^ + X^) (19) 

The factors p and g"^ in Eqs. (18,19) are defined in (12). In Fig. 7 it is verified that Eq. (10) with C = Coo indeed 
describes the numerical EG0E(l-|-2) results. Having demonstrated this, it is easily seen that the obvious condition 
for 5*'"*^° and PR (also strength functions) to be same in both h{l) and 1^(2) basis is 

Co(Ad) - Coo(Ad) =^ Ad = |A//.g| , C'(Ad) - 0.5 (20) 

Using Fig. 3 and the condition C^(A(i) ~ 0.5 gives for the m = 6, = 12 example, A^ = 0.29. In the dilute limit, the 
m dependence of A^ follows from Eqs. (12,13,20), 

Xd ^ A/(3to)1/2 (21) 

Thus the Gaussian domain arguments give v (in Xd ^ \/my) to be 0.5 unlike the improved BW domain arguments 
(Eq. (17)) giving 0.375. With Ad defined, a much more significant result that follows from (18-20) is Coo (A) = Co(Ad/A) 
and thus there is a duality in EG0E(l-|-2), i.e. the results in h{l) and V{2) basis are related to each other by the 
duality transformation A A^/A. As stated before, the same transformation is also derived by JV but using and 
^2 in the BW domain and this points-out its general validity. Strictly speaking Ad does not lie in the BW domain 
nor deep into the Gaussian domain. The duality transformation is well tested in Fig. 8 for Fk{E) and in Fig. 9 
for S™^°{E). In these calculations Ad = 0.29. It should be recognized that for the strength functions in Fig. 8, 
the variances are Co s-nd Cm in the h{l) and V{2) basis respectively. In the case of 5'™^°(£') one sees (from Fig. 9) 
departures, for Fk{E) close to Gaussian, in the region well away from the centroid of E and this could be because 
the tails of Fk{E) display exponential localization jM, 33] . These disagreements are not seen in as in this work 
only S™^°{E — 0) and £,2{E = 0) are studied. It is useful to point-out that there appears to be a close relationship 
between Ad and thermodynamics of finite quantum systems. Using the Gaussian domain formulas (see [20|') for the 
thermodynamic, information and single particle entropies, it is easily verified that at and around Ad, all the three 
entropies will be very close to each other; numerical verification of this result is given in po| . Therefore it is possible 
to define the region around Ad as 'thermodynamic region' for interac ting particle systems as here different definitions 
of thermodynamic quantities like entropy will give same results; see [23, . 
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VI. CONCLUSIONS 



In this paper an attempt is made to bring completion to the analytical (in BW and Gaussian domains) and 
numerical investigations, initiated by a number of research groups, of EG0E(l+2) random matrix model for finite 
interacting quantum systems. Towards this end, a function describing the BW to Gaussian transition in strength 
functions is identified (Eq. (4)) and it is used to study participation ratio and information and structural entropy 
as a function of the interaction strength. Also it is shown, using Gaussian domain results, that the duality point Xd 
behaves more like ^ where m is number of fermions. Applications of these results are given for the BW to 

Gaussian transition in the series of neutral atoms Cel, PrI, Ndl, Pml and SmI. As for EG0E(l+2), what remains is 
a rigorous analytical treatment of this random matrix model. This will give for example a theory for a vs A (see Fig. 
2), a theory for 5*'"^° and PR in the A ^ Ac domain etc. Finally it is important to remind that only recently rigorous 
analytical treatment has started becoming available for the simpler EG0E(2) [23l |. 
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FIG. 1: Strength functions Fk{E) for a 20 member EG0E(H-2) for various values of the interaction strength A in {H} — 
h{l) + X{V{2)} for a system of 6 fermions in 12 single particle states; the matrix dimension is 924. The single particle energies 
used in the calculations are — {i+ = 1, 2, . . . , 12 just as in 1]. In the figures Fk{E) is plotted against E — {E — e)/cr 

where e is the spectrum centroid and a is the width. The histograms are EG0E(l+2) results and the continuous curves are the 
best fit Fk:BW-g{E) from Eq. (4). In constructing the strength functions, |C^p are summed over the basis states |fe) in the 
energy window ± A and then the ensemble averaged Fj^^ (E) vs E is constructed as a histogram; the value of A is chosen 

to be 0.05 for A = 0.06 and beyond this A = 0.1. Here Ek ~ {Ek — ^h)/o-h and in the figures results shown for Fk{E) with 
Ek ~ 0. Similar results are also obtained for the (m — 7,N = 14) system. 
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FIG. 2: For A > 0.08, a vs A obtained by fitting the Fk{E) in h{l) basis to the interpolating form Fk;Bw-g given by (4). 
Results arc shown for the EG0E(l+2) system used in Fig. 1. The filled circles give the best fit a values and the continuous 
curve, given by a = 24/[exp — (40(A — Ao)) + 1] + 2 with Ao = 0.29, guides the eye. It is curious to note that Ao is close to Ad, 
the duality pont discussed in Section 5. Similar results are also obtained for the (m = 7, TV = 14) system. 




FIG. 3: Square of the correlation coefficient Co vs A for the EG0E(l+2) system used in Fig. 1. 
curve) given by (11,12) are compared with the numerical EGOE(l+2) results (open circles). 
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FIG. 4: Participation ratio 52(-E = 0)7^2 vs A for the EG0E(l+2) system used in Fig. 1. Theoretical results given by (14) 
(open squares) are compared with the EGOE(l+2) results (filled circles). For comparison, the Gaussian domain result (dashed 
curve) from (10) is also shown. 
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FIG. 5: Structural entropy Seti{E = 0) vs A for the EGOE(H-2) system used in Fig. 1. Theoretical results given by (14,8) 
(open squares) are compared with the EGOE(l+2) results (filled circles). For comparison, the Gaussian domain result (dashed 
curve) from (16) is also shown. 
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FIG. 6: Strength functions Fk{E) for Ccl to SmI. Histograms are calculated strength functions and the smooth curves are 
Fk:BW-g{E). Also given in the figure are the calculated 71 (skewness) and 72 (excess) values and the deduced values, from the 
best fits, of a characterizing Fk:Bw-g{E). 
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FIG. 7: exp(S''°''°(£; = 0) - S'g'ois) in V(2) basis (filled circles) for four different A values for the EGOE(l+2) system used in 
Fig. 1. The continuous curve is from Eq. (10) with ^ = ^^o- 
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FIG. 8: Strength functions Fk{E) in the and V{2) basis for four A values related by the duality transformation A A^/A. 
Results are for the EG0E(l+2) system used in Fig. 1. Here \d = 0.29. Similar results for the BW spreading widths are given 
in [l9| for several EGOE(l+2) systems with h{l) also chosen to be random. 
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FIG. 9: Same as Fig. 7 but for exp(S'"'°(S) 
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Similar results for ^2 but only at iJ = are given in 



